Spin-orbit coupling in curved graphene, fullerenes, nanotubes, and nanotube caps. 



Daniel Huertas-Hernando,^ F. Guinea,^ and Arne Brataas^^^ 

^Department of Physics, Norwegian University of Science and Technology, N-7491, Trondheim, Norway 
^Instituto de Ciencia de Materiales de Madrid, CSIC, Cantoblanco E28049 Madrid, Spain 
'^Centre for Advanced Study, Drammensveien 78, 0271 Oslo, Norway 

A continuum model for the effective spin orbit interaction in graphene is derived from a tight- 
binding model which includes the n and a bands. We analyze the combined effects of the intraatomic 
spin orbit coupling, curvature, and applied electric field, using perturbation theory. We recover the 
effective spin-orbit Hamiltonian derived recently from group theoretical arguments by Kane and 
Mele. We find, for fiat graphene, that the intrinsic spin-orbit coupling Aint cx: and the Rashba 
coupling due to an perpendicular electric field £, Ag oc A, where A is the intraatomic spin-orbit 
coupling constant for carbon. Moreover we show that local curvature of the graphene sheet induces 
an extra spin-orbit coupling term Acurv oc A. For the values of £ and curvature profile reported 
in actual samples of graphene, we find that Aint < Ae < Acurv. The effect of spin orbit coupling 
on derived materials of graphene like fullerenes, nanotubes, and nanotube caps, is also studied. For 
fullerenes, only Aint is important. Both for nanotubes and nanotube caps Acurv is in the order 
of a few Kelvins. We reproduce the known appearance of a gap and spin-splitting in the energy 
spectrum of nanotubes due to the spin-orbit coupling. For nanotube caps, spin-orbit coupling 
causes spin-splitting of the localized states at the cap, which could allow spin-dependent field-effect 
emission. 



INTRODUCTION. 



A single layer of carbon atoms in a honeycomb lattice, graphene, is an interesting two-dimensional system due to 
its remarkable low energy electronic properties^, |3, Q, £-9- a zero density of states at the Fermi level without an 
energy gap, and a linear, rather than parabolic, energy dispersion around the Fermi level. The electronic properties of 
the many realizations of the honeycomb lattice of carbon such, e.g. bulk graphite (3D), carbon nanotube wires (ID), 
carbon nanotube quantum dots (OD), and curved surfaces such as fullerenes, have been studied extensively during the 
last decade. However, its two dimensional (2D) version, graphene, a stable atomic layer of carbon atoms, remained for 
long ellusive among the known crystalline structures of carbon. Only recently, the experimental realization of stable, 
highly crystaUine, single layer samples of graphene 0, IMS Hi have been possible. Such experimental developments 
have generated a renewed interest in the field of two dimensional mesoscopic systems. The peculiar electronic properties 
of graphene are quite different from that of 2D semiconducting heterostructures samples. It has been found that the 
integer Hall effect in graphene is different than the "usual" Quantum Hall effect in semiconducting structures [sl Isl 
id, 11]. Moreover, it has been theoretically suggested that a variety of properties, e.g. weak (anti)localisation I4 
13 . l^^l^, shot-noise [1^ and anomalous tunneling-Klein's paradox^lgl, are qualitatively different from the 

behavior found in other 2D systems during the last decades. All these predictions can now be directly investigated 
by experiments. The activity in graphene, both theoretically and experimentally, is at present very intense. However 
so far, the work has mainly focused on i) the fact that the unit cell is described by two inequivalent triangular 
sublattices A and B intercalated, and ii) there are two independent fc-points, K and K' , corresponding to the two 
inequivalent corners of the Brillouin zone of graphene. The Fermi level is located at these K and K' points and 
crosses the tt bands of graphene (see Fig.pP for details). These two features provide an exotic fourfold degeneracy of 
the low energy (spin-degenerate) states of graphene. These states can be described by two sets of two-dimensional 
chiral spinors which follow the massless Dirac-Weyl equation and describe the electronic states of the system near 
the K and K' points where the Fermi level is located. Neutral graphene has one electron per carbon atom in the tt 
band, so the band below the Fermi level is full (electron- like states) and the band above it is empty (hole- like states). 
Electrons and holes in graphene behave like relativistic Dirac fermions. The Fermi level can be moved by a gate voltage 
underneath the graphene samplej^J. State-of-the-art samples are very clean, with mobilities fi ~ 15000cm^y~ls~l[6j, 
so charge transport can be ballistic for long distances across the sample. From the mobilities of the actual samples, 
it is believed that impurity scattering is weak. Furthermore, it has been recently suggested that the chiral nature of 
graphene carriers makes disordered regions trans par ent for these carriers independently of the disorder, as long as it 
is smooth on the scale of the lattice constant ll^l Il4l \2U^ . 

Less attention has been given to the spin so far. The main interactions that could affect the spin degree of 
freedom in graphene seem to be the spin-orbit coupling and exchange interaction. It is not known to which extent 
magnetic impurities are present in actual graphene samples. Their effect seem small though, as noticed recently 
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when investigating weak localization and universal conductance fluctuations in graphene|l4j. Spin-orbit interaction in 
graphene is supposed to be weak, due to the low atomic number Z = 6 of carbon. Therefore both spin splitting and 
spin-flip due to the combination of spin-orbit and scattering due to disorder is supposed to be not very important. 
As a result, the spin degree of freedom is assumed to have a minor importance and spin degenerate states are 
assumed. Besides, the spin degeneracy is considered to be "trivial" in comparison to the fourfold degeneracy previously 
mentioned, described by a pseudo-spin degree of freedom. At present , th ere is a large activity in the study of the 
dynamics of this pseudo-spin degree of freedomH |l[ll[Il[ll[ll[Illll[ll[ll[lllil|23IIl 

We think that the physics of the electronic spin in graphene must be investigated in some detail, however. Although 
it could be that the electronic spin is not as important/exotic as the pseudo-spin when studying bulk properties, edge 
states may be quite different. Induced magnetism at the edges of the surface of graphite samples irradiated with 
protons have been reported 27] . Moreover, perspectives of spintronic applications in graphene could be very promising, 
so it is important to clarify the role of the electronic spin. This is one of the main purposes of the present paper. 
Moreover, we feel that the existent knowledge about the spin-orbit interaction in graphene is not yet complete[2J| and 
that certain, both quantitative and qualitative, points must be discussed in more detail. That is why we focus our 
discussion on the spin-orbit coupling. The effect of other interactions as the exchange interaction will be discussed 
elsewhere. 

Spin-orbit coupling in graphene has an intrinsic part, completely determined from the symmetry properties of the 
honeycomb lattice. This is similar to the Dresselhauss spin-orbit interaction in semiconducting heterostructuresfis) . 
Group theoretical arguments allow to obtain the form of the effective Hamiltonian for the intrinsic spin-orbit coupling 
around the K, K' points I^ES 113 ■ predicted that this interaction opens-up a gap Aint in the energy dispersion. 

However, the strength of this intrinsic spin-orbit coupling is still a subject of discussion, although it is believed to be 
rather small, due to the weakness of the atomic intraatomic spin orbit coupling of carbon A. If an electric field £ 
is applied perpendicular to the sample, a Rashba interaction 31] A^ will also be present in graphene. Analogous to 
the intrinsic coupling, group theoretical arguments allow deducing the form of the Rashba interaction 24, 25]. The 
strength of this Rashba spin-orbit coupling is also still under discussion. 

We follow a different approach. We set up a tight binding model where we consider both the tt and a bands of 
graphene and the intraatomic spin-orbit coupling A. We also include curvature effects of the graphene surface and 
the presence of a perpendicular electric field £. Starting from this model, we obtain an effective Hamiltonian for 
the TT bands, by second order perturbation theory, which is formally the same as the effective Hamiltonian obtained 
previously from group theoretical met hods [2^ Isil] by Kane and Meleji^]. Moreover, we show that curvature effects 
between nearest neighbor atoms introduce an extra term Acurv into the effective spin-orbit interaction of graphene, 
similar to the Rashba interaction due to the electric field Ag . We obtain explicit expressions for these three couplings 
in terms of band structure parameters. Analytical expressions and numerical estimates are given in Table p. 

We find that the intrinsic interaction Ajnt ~ lOmK is two orders of magnitude smaller than what was recently 
estimated (241]. Similar estimates for Aint have been reported recently ,22„ ^2, 33]. Moreover, we find that for typical 
values of the electric field as e.g. used by Kane and Melejl^] Af ~ 70mK. Similar discussion for Ag has appeared also 
recently . So spin-orbit coupling for flat graphene is rather weak. Graphene samples seem to have an undulating 
surface however |14j. Our estimate for the typical observed ripples indicates that Acurv ~ 0.2K. It seems that curvature 
effects on the scale of the distance between neighbouring atoms could increase the strength of the spin-orbit coupling 
at least one order of magnitude with respect that obtained for a flat surface. More importantly, this type of "intrinsic" 
coupling will be present in graphene as long as its surface is corrugated even if £ =0 when Af =0. 

The paper is organized as follows: The next section presents a tight binding hamiltonian for the band structure 
and the intraatomic spin-orbit coupling, curvature effects and a perpendicular electric fleld. Then, the three effective 
spin-orbit couplings Aint , Ag , Acurv for a continuum model of the spin-orbit interaction for the tt bands in graphene at 
the K and K' points are derived. Estimates of the values are given at the end of the section. The next section applies 
the effective spin-orbit hamiltonian to: i) fuUerenes, where it is shown that spin-orbit coupling effects play a small role 
at low energies ii) nanotubes, where known results are recovered and iii) nanotubes capped by semispherical fuUerenes, 
where it is shown that the spin orbit coupling can lead to localized states at the edges of the bulk subbands. The last 
calculation includes also a continuum model for the electronic structure of nanotube caps, which, to our knowledge, 
has not been discussed previously. A section with the main conclusions completes the paper. 
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DERIVATION OF CONTINUUM MODELS FROM INTRA ATOMIC INTERACTIONS. 

Electronic bands. 

The orbitals corresponding to the cr bands of graphene are made by Hnear combinations of the 2s, 2px and 2py 
atomic orbitals, whereas the orbitals of the vr bands are just the pz orbitals. We consider the following Hamiltonian: 

^ = ^SO + ^atom + ^TT + ^o- , (1) 

where the atomic hamiltonian in the absence of spin orbit coupling is: 

^atom — ^ ^ C^g/Cis' ^ ^ ^s^g'Cs^s^ • (2) 

i=x,y,z;s' = '\ ,1 s;s' = '\,l 

where ep.s denote the atomic energy for the 2p and 2s atomic orbitals of carbon, the operators Ci-s' and Cs-s' refer to Pz, 
Px, Py and s atomic orbitals respectively and s' denote the electronic spin. Tiso refers to the atomic spin-orbit 

coupling occuring at the carbon atoms and the terms TL-^^Tia describe the tt and a bands. In the following, we will 
set our origin of energies such that = 0. We use a nearest neighbor hopping model between the pz orbitals for 7i^, 
using one parameter T^p^ . The rest of the intraatomic hoppings are the nearest neighbor interactions Vppo- , Vgpa and 
Vssa between the atomic orbitals s,Px,Py of the a band. We describe the a bands using a variation of an analytical 
model used for three dimensional semiconductors with the diamond structure |34|. and which was generalized to the 
related problem of the calculation of the acoustical modes of graphene [l^- The model for the sigma bands is described 
in Appendix A. The bands can be calculated analytically as function of the parameters: 



V2 = ^ . (3) 

The band structure for graphene is shown in Fig.Q. 

Intraatomic spin-orbit coupling. 

The intraatomic spin orbit coupling is given by Tiso — ALs 36] where L and s are, the total atomic angular 
momentum operator and total electronic spin operator respectively, and A is the intraatomic spin-orbit coupling 
constant. We define: 
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FIG. 1: "(Color online)": Black (Full) curves: a bands. Red (Dashed) curves: vr bands. The dark and light green (grey) arrows 
give contributions to the up and down spins at the A sublattice respectively. The opposite contributions can be defined for the 
B sublattice. These interband transitions are equivalent to the processes depicted in Fig.^]. 



Using these definitions, the intraatomic spin-orbit hamiltonian becomes: 
The Hamiltonian Eq. © can be written in second quantization language as 



(5) 



(6) 



where the operators cl^y.^, and Cz^x,y;s' refer to the corresponding pz, Px and py atomic orbitals. The intraatomic 
hamiltonian is a 6 x 6 matrix which can be spht into two 3x3 submatrices: 

^so = 



nil 



(7) 



The block Tigo ^'^^^ ^^'^ basis states \pz t)i \Px i) and \py |): 



On the other hand Hgo- 



^so 



-1/22 




SO 




(8) 



(9) 
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FIG. 2: "(Color online)": Sketch of the relevant orbitals, and pz needed for the analysis of spin-orbit effects in a curved 
nanotube. The arrows stand for the different hoppings described in the text. 



acts on \pz i), \px T) and \py t) states. The eigenvalues of these 3x3 matrices are +A (J = 3/2) which is singly 
degenerate and —A/2 (J = 1/2) which is doubly degenerate. 

The term L+s_ + i_s+ of the intraatomic spin orbit coupling Hamiltonian, Eq.©, allows for transitions between 
states of the tt band near the K and K' points of the Brillouin zone, with states from the a bands at the same points. 
These transitions imply a change of the electronic degree of freedom, i.e. a "spin- flip" process. 

We describe the a bands by the analytical tight binding model presented in Appendix A. The six a states at 
the K and K' points can be split into two Dirac doublets, which disperse linearly, starting at energies ea-{K,K') — 
Vi/2 ± ^%V^I^ + and two flat bands at tcr{K, K') = -Vi ± V2. We denote the two a Dirac spinors as ipai and 
■00-2, and the two other "flat" orbitals as (pai and 00-2. The intraatomic spin orbit hamiltonian for the K and K' point 
becomes: 



T^SOK' = 



A 
2" 

sin 



{^°^ (?) 



*WT(f)V'<TiSKi(f) 



(I) 



'a 



*WT(^)^-2AKi(f)+«'^3;^^(r)0.2SKi(f)]} , yg 

^ {"^^^ (?) '^AK"ii^)i^'y^AK'i{r) + *^^,^(r)?/'^2SK'i(f) 

^ , ^ ^ 1 ~1 /2" r 

BK'l[r)\ I 



*Vt(i') + *W't(?) </'2i(?) + Mio) 



where stands for the two component spinor of the tt band, and cos(q!/2) and sin(Q;/2) are matrix elements given 
by: 



a — arctan 



(3^1 )/2 



(11) 



Next we would like to consider two posibilities: i) A curved graphene surface, ii) The effect of a perpendicular 
electric field applied to flat graphene. In the latter case we will have to consider another intraatomic process besides 
the intraatomic spin-orbit coupling, the atomic Stark effect. 
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Effects of curvature. 

In a curved graphene sheet, a hopping between the orbitals in the tt and a bands is induced (37j. First we wiU use 
a simple geometry to iUustrate the effect of curvature between neighbouring atoms. This geometry is schematically 
shown in Fig.|2], where we first consider two atoms at the same height along the axis of the tube. In this geometry 
we consider that the pz orbitals are oriented normal to the surface of the nanotube, the Px orbitals are oriented along 
the surface circumference (Fig.j2] ) and the Py orbitals are parallel to the nanotube axes. The curvature modifies the 
hopping between the two atoms compared to the flat surface for the pz and px orbitals but will not change, for this 
simple case, the hopping between py orbitals. The (reduced) Pz-Px hopping hamiltonian is the sum of three terms: 

Ht = ['^PP^ cos^{0) + Vppa sin2(6l)] 4i,,c^os' - [Vppyr sin2(6l) + Vppa cos^{9)] 4i,,c^os' + 

+ Vspa sin^(6')4^^,CsOs' + sin(6') cos(6') {Vpp-^ - Vppa) (clj^^^CxOs' ~ c^is-c^os') + h.c. (12) 

where and 1 denote the two atoms considered and 9 is the angle between the fixed Z axis and the direction normal 
to the curved surface (Fig.j2] ). The angle 9, in the limit when the radius of curvature is much longer than the 
interatomic spacing, a <C i?, is given by 6* « a/R. 

The hopping terms induced by (intrinsic) curvature discussed here break the isotropy of the lattice and lead to an 
effective anisotropic coupling between the tt and a bands in momentum space. 

The previous discussion can be extended to the case of general curvature when the graphene sheet has two different 
curvature radii, Ri and i?2 corresponding to the x and y directions in the plane. In that case, the factor R~^ has to 
be replaced by + i?^^. We now expand on ~ a/Ri,2 <C 1. By projecting onto the Bloch wavefunctions of the 
TT and a bands at the K and K' points, we find: 

Htk = (^pp- " ^pp-) + / '^'^{'^°Kf) [*AKT(^)^-is^T(f) + «'WT(i')^-i^^T(f) 

+ sin (I) \^\k^ (f)V'<T2SKT (f) + *LxT (f") V'^2AKT (f)] ]+h.c. (13) 

and a similar expression for the K' point. 

The induced spin orbit coupling, however, includes only contributions from the four a bands at K and K' with 
ta = ^i/2 ± \/ {9Vi)/A + V2 , as those are the only bands coupled to the tt band by the intraatomic spin orbit term 
considered here, Ea. (|10|l . We now assume that the energies of the a bands are well separated from the energy of the 
TT bands [e.^ = at the K and K' points). Then, we can use second order perturbation theory and from Ea. (|l()|l and 
Ea. p3|l we obtain an effective hamiltonian acting on the states of the tt band. 

Effect of an electric field. 

Now we discuss the atomic Stark effect due to a perpendicular electric field £. In this case, we need to consider the 
\s) orbital of the a bands at each site, and the associated hopping terms. The hamiltonian for this case includes the 
couplings: 

T-Lg ^ ^ (^Xe£c\^.^,Ciz-s' + ^sc\g.s''^is-s' +h.c'}j +Vspa ^ {axc\x.s''^Os-s' + ayc\y.,.,CQs-s' + h.c^ (15) 
i=l,2;s' = T,i s' = T,i 

where A = (pzlzjs) is a electric dipole transition which induces hybridization between the s and pz orbitals and 
where ax and ay are the x and y components of the vector connecting the carbon atoms and 1. First, we consider 
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the situation a^: = 1 and ay = 0. Again Vspa- is the hopping integral between the 2s and 2px,2py atomic orbitals 
corresponding to the cr band. We can now have processes such as: 

e 



(16) 

The intermediate orbitals |sO) and \px^) are part of the sigma bands. As before, we describe them using the analytical 
fitting discussed in Appendix A. The |sO) is part of the dispersive bands, and it has zero overlap with the two non 
dispersive a bands. The processes induced by the electric field, in momentum space, lead finally to: 



cos 



(f 



ba2BKl{f) I + h.C. 



(f) 



+ 



(17) 



Note that this hamiltonian mixes the states in the n band with states in the a bands which are orthogonal to 
those in Ec ljlOl) . Combining Eq.lQTf and Ea. (|10() we obtain, again by second order perturbation theory, an effective 
hamiltonian for the tt band: 



.2^2 ^Xe£V2 



2V2 AXe8V2 



3 2V, 



(18) 



The zero overlap between the states in the a bands in Ea. p7|) and Ea. (|10|l imply that only transitions between 
different sublattices are allowed. 
Defining a 4 x 4 spinor 



KiK') 



/*AT(f)\ 
*Ai(r) 

«'B|(f) 
V*Bi(f)/ 



(19) 



K{K') 



it is possible to join Eas. (|14|) and (|18|l in the following compact way: 



(20) 



where 



Ti-RK'n — 



1* 



K' 



As 



i£ + Ac 
AV2 



Ar. 



2Vl 
AVi 



2Vl 



2V2 



2V2 A\e£ 



^ppiz ) 



{-- 









Ar 

2 



1* 



K' 



a 



a 



(21) 



(22) 



where the limit Vi <C V2 (widely separated a bands) has been considered to approximate the above expressions. 

Eas. (|20l21l22|l constitute one of the most important results of the paper. First, we recover the effective form for the 
"Rashba-type" interaction expected from group-theoretical arguments recently |2J|. Even more importantly, our result 
shows that this effective spin-orbit coupling for the tt bands in graphene to first order in the intraatomic spin-orbit 
interaction A is given by two terms: 
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FIG. 3: "(Color online)": Sketch of the processes leading to an effective intrinsic term in the vr band of graphene. Transitions 
drawn in red (dark grey), and indicated by SO, are mediated by the intraatomic spin-orbit coupling. 



• Ag: Corresponds to processes due to the intraatomic spin-orbit coupling and the intraatomic Stark effect 
between different orbitals of the tt and a bands, together with hopping between neighboring atoms. The mixing 
between tlie tt and a orbitals occurs between the Pz and s atomic orbitals due to the Stark effect A and between 
the "Pz and Px.y due to the atomic spin-orbit coupling A. This contribution is the equivalent, for graphene, to 
the known Rashba spin-orbit interactionlsTj and it vanishes at E =0. 

• Acuiv: Corresponds to processes due to the intraatomic spin-orbit coupling and the local curvature of the 
graphene surface which couples the tt and a bands, together with hopping between neighboring atoms. The 
mixing between the tt and a orbitals in this case occurs between pz and "Px.y atomic orbitals both due to the 
atomic spin-orbit coupling A and due to the curvature. This process is very sensitive to deformations of the 
lattice along the bond direction between the different atoms where the p part of the the s]? orbitals is important. 



Intrinsic spin orbit coupling. 

We can extend the previous analysis to second order in the intraatomic spin-orbit interaction A. We obtain effective 
couplings between electrons with parallel spin. The coupling between first nearest neighbors can be written as: 

b.o T) ^ b.o i) ^ i) ^ \v^ T) 

b.Ot) ^ -ib.Oi)-H^Mi) ^ i|p^2i)-^|p,2i) ^Mt) (23) 

where the label stands for the central atom. These three couplings are equal, and give a vanishing contribution at the 
K and K' points. The intrinsic spin-orbit coupling vanishes for hopping between neighbouring atoms, in agreement 
with group theoretical arg umentsHliill^. We must therefore go to the next order in the hopping. Expanding 
to next nearest neighbors, we find a finite contribution to the intrinsic spin-orbit coupling in a flat graphene sheet, 
corresponding to processes shown in Fig.P]. 

In this case, both the dispersive and non dispersive bands contribute to the effective tt — tt coupling, as schematically 
shown by the different arrows in Fig.^. In order to estimate quantatively the magnitude of the intrinsic coupling, 
we consider processes represented in Fig.p], which are second order in A, in momentum space, finally obtaining: 

^ ,3A^ 

nintKiK') ^^y^ _ v^){2Vl + Vi) 

'^'f*L(K')T(^)*^^(^')T(?) - '^\K{K')ii'^)'^AK(K')l{?) - *1j,^(^,)i.(?)*bk(K')t(?) + ^BK(K')l{?)^BK(K')i{ 

(2 

where the ± sign corresponds to K(K') respectively. We define the intrinsic spin-orbit coupling parameter Aint in 
the limit Vi <C V2 (widely separated a bands) as: 



3 A2 Yl 3 A2 I'Vi 

4 W (Vi ~ v^){2v^ + Vi) ~ 1~v yvz 



^i"t- 7 — 7772 772777772-T-F727 - 7— ( 7F ) (25) 
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Intrinsic coupling: Aint 


3A^ (V^X 
4 \\ \V^j 


O.OIK 


Rashba coupling (electric field £ « 50V/300nm): 


2V2 AAe£ 
3 V, 


0.07K 


Curvature coupling: Acurv 


A{Vpp„-Vpp^) (a 1 a.\ (ViY 
Vi \Rl ^ R2 I \ V2 I 


0.2K 



TABLE I: Dependence on band structure parameters, curvature, and electric field of the spin orbit couplings discussed in the 
text in the limit Vi <C V2 (widely separated a bands). The parameters used are A « 0.264A36], £ ^ SOV/SOOnm"?, , 

A = 12meV [H, 39],Kp<T ~ 4.2eV, Vss<t 3.63eV, Vpp^ ~ 5.38eV and Vpp-^ 2.24eV |40„i41J, Vi = 2.47eV, V2 = 6.33eV, 

a = I.42A, Z ~ 100 A, ft ~ 10 A and i? ~ 50 - lOOnm. 



Our Hamiltonian, Ea. H24|) . is equivalent to the one derived in p4 



HtSO intrinsic = j <f'?^\nt'^^[TzCrzS,]^ 



(26) 



where fz — ±1 denotes the K{K') Dirac point and ^ — (^k, k'Y 



NUMERICAL ESTIMATES. 



We must now estimate Aint, and Acurv We have A = ?>ao/Z « 0.264A|3^, where Z = 6 for carbon and a„ i s 
the Bohr radius, £ ^ SOV/SOOnmjliil, the atomic spin-orbit splitting for carbon A = 12meV 1.3 x lO^KjUlHil, 
the energy difference between the 7r-2pz orbitals and the cr-sp^ orbitals e^^ — €a ^ (14.26 — 11.79)eV= 2.47eV, the 
energy difference between the 2p and the 2s atomic orbitals eg — ep ^ (19.20 — 11.79)eV= 7.41eV and the hoppings 



between the 2s, 2 
eV - 



2pz orbitals of neighbouring atoms as Vspa- ~ 4.2eV, Vgscr — 3.63eV and Vppa ~ 5.38eV and 
2.24eV |40|,|41|. We have Vi = 2.47eV and V2 = 6.33eV. We estimate Ai„t ~ {3A^/4yi){Vi/V2)'^ ~ 0.1 x lO'^ 
O.OIK. Aint is two orders of magnitude smaller than the estimate in[2^. The discrepancy seems to arise first, 
because the intrinsic spin-orbit splitting Aint estimated here is proportional to the square of the intraatomic spin-orbit 
couphng A^, instead of being proportional to it, as roughly estimated inQ. Besides, a detailed description of the a 
bands is necessary to obtain the correct estimate. Spin-orbit splittings of order 1-2.5 K have been discussed in the 
literature for graphite |42l|. However in graphite, the coupling between layers is important and influences the effective 
value of the spin-orbit splitting, typically being enhanced with respect to the single layer value [30I li^ lil] . 

For the other two couplings, we use the full expression obtained for Af and Acurv and not the limiting form 
< 1^2, in order to be as accurate as possible. First, we obtain As = (2\/2/3) {\e£AV2/ (2Vf + V^)) ^ 0.6 x 10^'' 
eV 0.07K. This estimate for Ag, depends on the external electric field chosen. Our estimate, for the same value 
of the electric field, is two orders of mag nitude bigger than the estimate inlH. So far curvature effects have been 
excluded. Curvature effects will increase the total value for the effective spin-orbit interaction in graphene. Graphene 
samples seem to have an undulating surface |l4l| . The ripples observed seem to be several A height and a few tens nm 
laterally 01 . First we consider the simplest example of a ripple being a half-sphere of radius R. The part of the sphere 
which intersects the plane of flat graphene and therefore constitutes the ripple, is assumed to have a typical height 
h ^ R so the radius R is roughly of the same order of ma gnit ude of the lateral size in this case. It seems pos sible 
to identify ripples of lateral size ranging 50nm -lOOnm in [lj|. Choosing a = 1.42A and Ri ^ R2 ^ 50nmll4j, we 
obtain Acurv = {2a/ R) {{Vpp„ - Vpp^) AVi/ (2Vf + Vi)) ~ 2.45 x lO-^eV ^ 0.28 K. Choosing Ri - R2 - lOOnm 
we obtain Acurv ~ 1-22 x 10~^eV 0.14 K. Now we consider a different model where we assume that the sample 
has random corrugations of height h and length I [l^ . The graphene surfacepresents then an undulating pattern of 
ripples of average radius R ~ f/h^. Choosing I ~ lOOA and /i ~ 10 AO, we obtain again R ^ lOOnm, which 
leads to the same value for Acurv ^ 1-22 x 10~^eV — > 0.14 K. In any case, it seems clear that due to curvature effects, 
the effective spin-orbit coupling in graphene could be higher for curved graphene than for perfectly flat graphene. 
Moreover, spin-orbit coupling in (curved) graphene would be present even for £ =0. A more detailed discussion/study 
of the local curvature/corrugation of graphene is needed in order to obtain more accurate estimates. 

To conclude this section we present the effective hamiltonian for the 7r-bands of graphene including the spin-orbit 
interaction: 



fz&ydy 



] + Ai 



int TzCTzSz 



A^ 
2 ■ 



J y Ox 



(27) 



10 



where h vF = •v{ 37na/2, a ~ 2.46 A being the lattice constant for graphene, 70 ~ 3eV the McClure intralayer couphng 
constant pui 1451 ■ Aint ~ O.QIK, Ar = Ag + Acuiv the Rashba-Curvature coupling (RCC), where Ag ~ Q.07K for £ w 
50V/300nm and Acurv ~ 0.2K. Table[l] summarizes the main results obtained in the paper for the effective spin-orbit 
couplings in a graphene layer. 



APPLICATION TO FULLERENES, NANOTUBES AND FULLERENE CAPS. 

Spherical fuUerenes. 

When topological deformations in the form of pentagons are introduced in the hexagonal lattice of graphene, curved 
structures form. If 12 pentagons are introduced, the graphene sheet will close itself into a sphere forming the well 
know fuUerene structure . The usual hexagonal lattice "lives" now on a sphere and presents topological defects in 
the form of pentagons. The continuum approximation to a spherical fuUerene leads to two decoupled Dirac equations 
in the presence of a fictitious monopole of charge ±3/2 in the center of the sphere which accounts for the presence of 
the 12 pentagons'??, 'is^. The states closest to the Fermi level are four triplets at e = 0. We consider the effect of the 
spin-orbit coupling on these triplets first. 

Both the coupling induced by the Rashba-Curvature, Ea. (|20l I^Tl) . and the intrinsic coupling, Ea. (|24ll . are written 
in a local basis of wavefunctions where the spin is oriented perpendicular to the graphene sheet, \9, 0, J-t)i l^j -Li)- 
It is useful to relate this local basis with a fixed basis independent of the curvature. Such relation depends on the 
curvature of the graphene sheet considered. In the case of a fuUerene the graphene sheet is in a sphere. The details 
are given in Appendix B. 

The gauge field associated to the presence of fivefold rings in the fuUerene can be diagonalized using the basis: 

^AKks(?) = *AXks(f) + ^*SA"k.(f) 

^BK'Ws *SK'ks (r ) + *^A'k. (f ) ■ (28) 

Equivalent transformation is obtained exchanging A B. 

In this basis, the wavefunctions of the zero energy states are|47): 

l-fls/C) = 

\OsJC) = 

|-ls/C) = 

l + ls/C') = 

|Os/C') = 

where \AK) and \BK') are envelope functions associated to the K and K' points of the Brillouin zone and corre- 
sponding to states located at the A and B sublattices respectively. Note that, at zero energy, states at K{K') are 
only located at sublattice A{B) sites. \s) denotes the usual spinor part of the wave function corresponding to the 
electronic spin s The hamiltonian Tiint couples orbitals in the same sublattice whereas TYr couples orbitals in 

different sublattices. So Tip; has zero matrix elements between zero energy states, as it does not induce intervalley 
scattering, mixing K and K' states |13j. 
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In the {| + 1 t), I + 1 J,), |0 t), |0 J,), I — 1 t), I — 1 i)} basis, the Hamiltonian for the K. point of a fuUerene looks Hke: 



/ Aint 





'T-S-Oint 




^int 

V2Aint 









y2Ai„t 



\ 










V 



The Hamiltonian for K,' is '^s-Oint — ~^s-Oinf 



V2Ai„t -Aint 

Aint ) 



(30) 



Diagonalizing the Hamiltonian Eq. ()30(l , we obtain that each set of spin degenerate triplets obtained in the absence 
of the spin-orbit interaction split into: 



e = +Aint- 
e = -2Aint 



*a:{| + 1T),|-U),V:7I + U) + 



^|OT),/^|-iT) + 



Oi)} 



*-2A:{7| + U)-ji|0T), 



^|-1T)-J>i)} 



(31) 



Each of these solutions is doubly degenerate, corresponding to the JC and K! points. In principle, many body effects 
associated to the electrostatic interaction can be included by following the calculation discussed inj49i|. 



Spin-orbit coupling in nanotubes. 

The previous continuum analysis can be extended to nanotubes. We use cylindrical coordinates, 0, (}>, and, as before, 
define the spin orientations | \ [) as parallel and antiparallel to the z axis. The matrix elements relevant for this 
geometry can be easily obtained from Ea. (|55|l in Appendix B, by choosing = -k 12. The eigenstates of the nanotube 
can be classified by longitudinal momentum, fc, and by their angular momentum n, e±^k,n = ±fiwF\/fc^ + n^/R^, where 
R is the radius of the nanotube. After integrationg over the circumference of the nanotube / dcj), the Hamiltonian of 
a nanotube including spin-orbit interaction is: 

/ \At) \_f hvF{k- in/R) + T iArtts, \ ( \At) \ 

ns-OR 1^^^ J - y j^^^^j^ ^ _ ^ iAflTTs, )\ \Bt) ) ^"^^^ 

where the r — ±1 corresponds to the K{K') Dirac point. Note the basis states \At) and \Bt) used to define Ea. H32|l 
are spinors in spin subspace where the matrix Sz acts on (see Appendix C for details). The contribution from the 
intrinsic spin-orbit Ajnt becomes zero after integrating over the nanotube circumference (Appendix C) . The spin-orbit 
term iA/jTrSj in Ea. (|32|l is equivalent to the term proportional to Uy obtained in Eq. (3.15, 3.16) of Ref. [37|. It is 
important to note that the spin orientations | \ [) in Ea. (|32|) are defined along the nanotube axis, whereas the spin 
orientations used in Eq. (3.15, 3.16) of Ref. (33| are defined perpendicular to the nanotube surface. On the other 
hand, we do not find any contribution similar to the term proportional to Ux{r) in Eq.(3.15, 3.16) in Ref. |37l |. In any 
case, such contributions are not important as they vanish after integrating over the circumference of the nanotube|37l|. 
Besides, our results are in agreement with the results obtained in jSOj. 

The energies near the Fermi level, n = 0, are changed by the spin-orbit coupling, and we obtain: 



ek = ±^{nA^f + {hvFkY. (33) 



There is an energy gap ttAr at low energies, in agreement with the results in[33,|51|. The ttAr, gap originates as a 
consequence of the Berry phase gained by the electron and hole quasiparticles after completing a closed trajectory 
around the circumference of the nanotube under the effect of spin-orbit interaction Ar | 37j . Similarly, Ar will give 
rise to a small spin splitting for ?i 7^ [s^, |^ 

Cfc = ±\j (^Ar)' + (?iWF)2(fc2 + {n/Rf) + 2{72/R)hvFARns,. (34) 

For a single wall nanotube of radius Ri ^6,12, 24 A and i?2 — > 00, a ^ I.42A and for £ =0, we get Ar 12, 6, 3K 
respectively. 
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a) 



b) 



FIG. 4: "(Color online)": Left: One fifth of a fuUerene cap closing an armchair nanotube. The full cap is obtained by gluing 
five triangles like the one in the figure together, forming a pyramid. A pentagon is formed at the apex of the pyramid, from 
the five triangles like the one shaded in green (grey) in the figure. The edges of the cap are given by the thick black line. The 
cap contains six pentagons and seventy hexagons, and it closes a 25 x 25 armchair nanotube. Right: Sketch of the folding 
procedure of a fiat honeycomb lattice needed to obtain an armchair nanotube capped by a semispherical fuUerene. 



Nanotube caps. 



Localized states at zero energy. 



An armchair (5A^ x bN) nanotube can be ended by a spherical fuUerene cap. The cap contains six pentagons and 
5N{N + l)/2 — 5 hexagons. When = 3 x M, the nanotube is metalhc. An example of such a fuUerene cap is given 
in Fig.^. The boundary between the semispherical fuUerene and the nanotube is a circle (Fig.[S]). The solutions of 
the continuum equations have to be continuous accross this boundary, and they have to satisfy the Dirac equations 
appropriate for the sphere in the cap and for the torus in the nanotube respectively. 

The boundary of the nanotube in the geometry shown in Fig. [3] is a zigzag edge. Hence, zero energy states can be 
defined [53, 113, which at this boundary will have a finite amplitude on one sublattice and zero on the other. There is 
a zero energy state, |^n) at this boundary, for each value of the angular momentum around the nanotube n. They 
decay towards the bulk of the nanotube as 

*„(z,0,if') = c'e^"0e-(n^)/fl. „<o (35) 

where we are assuming that the nanotube is in the half space z > (see Fig.jS]). 

A zero energy state in the whole system can be defined if there are states inside the gap of the nanotube ttAr, which 
can be matched to the states defined in Ea. H35|) . At the boundary we have 9 = 7r/2, cos(6'/2) = sin(6l/2) = l/\/2. 
Hence, we can combine states \l s JC) and \l s K.'), I — ±1 in Ea. (|29() in such a way that the amplitude at the boundary 
on a given sublattice vanishes: 



V2 



(| + ls/C) + | + ls/C')) 



Is) 



V2 



(|-ls/C) + |-ls/C')) 



+ 1.s)b^^(| + 1s/C)-| + 1./C')) = 



8n 



\AK) 


\AK) 






i\BK') 


i\BK') 



(36) 



These combinations match the states decaying into the nanotube, Ea. l|35|l . This fixes the constant C in Eq.JSHJ to 
be C = ■\/3/(87r). Note that the wavefunctions with I = can only be matched to states that not decay into the bulk 
of the nanotube, i.e. with n = 0. 
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7. 



FIG. 5: "(Color online)": (Up) Sketch of the matching scheme used to build a zero energy state at a fuUerene cap. (Down)The 
wavefunction is one half of a zero energy state at the cap, matched to a decaying state towards the bulk of the nanotube. See 
text for details. 



Thus, there are two states per spin s, | + 1 s), | — 1 s), localized at the cap and with finite chirality, n = ±1. A 
sketch of this procedure is shown in Fig.jS]. This continuum approximation is in general agreement with the results 
in_54]. 

As in the case of a spherical fuUerene, only the intrinsic spin-orbit coupling mixes these states. The energies of the 
I ± 1 s)y!i states are not affected by Ajnt, as the contribution from | + 1 s), is cancelled by the contribution from | — 1 s). 
On the other hand, the states | ± 1 s)s split in energy, as the | + 1 s) and | — 1 s) contributions add-up: 

| + l,s=t,i>B^eTa =±Aint 

|-l,s=r,i)s->eT,i = TAi„t. (37) 
Note that each state has a finite chirality. 



Localized states induced by the spin-orbit interaction. 

The remaining states of a spherical fuUerenes have multiplicity 21 + 1, where I > 2 and energy e; = 
■iihv-p / R\J l{l + 1) — 2. The angular momentum of these states along a given axis, m, is — Z < m < /. The sub- 
bands of the nanotube with angular momentum ±m, have gaps within the energy interval — A„i = —hvp\fn\/ R < e < 
Am = v-p\m\/R. Thus, there is a fuUerene eigenstate with I = 2 and angular momentum m = ±2 which lies at the 
gap edge of the nanotube subbands with the same momentum. The fuUerene state is: 

^ (n ^ ( sin(0) [1 + cos(0)] \K) - [1 + cos{e)f\K') \ 

^i=ira=2[pA)- l^i[n_cos(0)]2|X)+isin(0)[l + cos(0)]|i^') y ^ ' 

which can be matched, at 6* = 7r/2, to the nanotube eigenstate: 

1 ^2.0^ \K)-\K') 



The spin orbit coupling acts as a position dependent potential on this state, and it shifts its energy into the m — 2 
subgap, leading to the formation of another localized state near the cap. 

In the following, we consider only the Rashba-Curvature coupling Aji. In order to analyze the extension of the 
state, we assume that the localized state decays in the nanotube, z > as: 

^ ^(z 6) = ^ ^2i<t> ^.z/R ( \K) - \K') \ 
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We match this wavefunction to that in Ea. H39|l muhiphed by the same normahzation constant, C. We assume that 
the state is weakly locahzed below the band edge, so that the main part of the wavefunction is in the nanotube, 
and K ^ 1. Then, we can neglect the change in the spinorial part of the wavefunction, and we will fix the relative 
components of the two spinors as in Ea. (|39|l and Ea. (|40|l . The normalization of the total wavefunction implies that: 

where the first term in the r.h.s. comes from the part of the wavefunction inside the cap, Ea. H39|) . and the second 
term is due to the weight of the wavefunction inside the nanotube, Ea. (|40|l . As expected, when the state becomes 
delocalized, k — s- 0, the main contribution to the normalization arises from the "bulk" part of the wavefunction. 
The value of n is fixed by the energy of the state: 

where the energy of the state e is inside the subgap of the nanotube because of the shift induced by the spin orbit 
interaction A^. 

We now calculate the contribution to the energy of this state from the Rashba-Curvature coupling, which is now 
finite, as this state has weight on the two sublattices: 

The r.h.s. in Ea. l43() can be described as the sum of a bulk term, ±Ar/4, and a term due to the presence of the 
cap, whose weight vanishes as the state becomes delocalized, k ^ 0. The absolute value of the Rashba-Curvature 
contribution is reduced with respect to the bulk energy shift, which implies that the interaction is weaker at the cap. 

The bulk nanotube bands are split into two spin subbands which are shifted in opposite directions. The surface 
states analyzed here are shifted by a smaller amount, so that the state associated to the subband whose gap increases 
does not overlap with the nanotube continuum. Combining the estimate of the energy between the gap edge and the 
surface state in Ea. (|43|l and the constraint for n in Ea. H42|) . we find: 

59ARi? 

Finally, the separation between the energy of the state and the subgap edge is: 

(45) 

For a Ceo fuUerene of radius R ^ 3.55A we obtain, for £ =0, a value Ai^/4 ^ 3K. This effect of the spin orbit 
interaction can be greatly enhanced in nanotube caps in an external electric field, such as those used for field emission 
devices In this case, spin-orbit interaction in may allow for spin dependent field emission of such devices. The 
applied field also modifies the one electron states, and a detailed analysis of this situation lies outside the scope of 
this paper. 



CONCLUSIONS. 



We have analyzed the spin orbit interaction in graphene and similar materials, like nanotubes and fuUerenes. We 
have extended previous approaches in order to describe the effect of the intraatomic spin orbit interaction on the 
conduction tt and valence a bands. Our scheme allows us to analyze, on the same footing, the effects of curvature 
and perpendicular applied electric field. Moreover, we are able to obtain realistic estimates for the intrinsic Ajnt 
and Rashba-Curvature /S.^. = Ag + Acurv effective spin-orbit couplings in graphene. We have shown that spin-orbit 
coupling for fiat graphene is rather weak Ajnt ~ lOmK and Ag ^ 70mK for f =50V/300nm. Moreover curvature 
at the scale of the distance between neighbouring atoms increases the value of the spin orbit coupling in graphene 
Acurv ^ Af ^ Aint. This is because local curvature mixes the tt and a bands. Graphene samples seem to have an 
undulating surface [lj|. Our estimate for the typical observed ripples indicates that Acurv could be of order ~ 0.2K. 
A more detailed study of the curvature of graphene samples is needed in order to obtain a more precise estimate. 
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We conclude that the spin-orbit couphng A^j expected from symmetry arguments |2^, has an curvature-intrinsic part 
besides the expected Rashba couphng due to an electric field = Ag -|- Acurv Therefore A^j can be higher for curved 
graphene than for flat graphene. To our knowledge, this is the first time that this type of "new" spin-orbit coupling 
has been noticed. A^urv is in a sense a new "intrinsic/topological" type of spin-orbit interaction in graphene which 
would be present even if £ = 0, as long as the samples present some type of corrugation. One important question 
now is how these ripples could affect macroscopic quantities. It has been already suggested that these ripples may 
be responsible for the lack of weak (anti)localization in graphene [l3lll3|. Other interesting macroscopic quantities 
involving not only the pseudo spin but also the electronic spin maybe worth investigation. These issues are beyond 
the scope of the present paper and will be subject of future work. 

It is also noteworthy that our estimates A^ ^ Ajnt, is opposite to the condition A^j ^ Ajnt obtained by Kane 
and Meleli^] to achieve the quantum spin Hall effect in graphene. So the quantum spin Hall effect may be achieved 
in neutral graphene {£ = 0) only below ^ O.OIK and provided the sample is also free of ripples so the curvature 
spin-orbit coupling Acurv ^ Ajnt- Further progress in sample preparation seems needed to achieve such conditions 
although some preliminary improvements have been recently reDorted|l4l|. Moreover, corrugations in graphene could 
be seen as topological disorder. It has been shown that the spin Hall effect survives even if the spin-orbit gap Ajnt is 
closed by disorder |2^. A detailed discussion of the effect of disorder on the other two spin-orbit couplings Ag, Acurv 
will be presented elsewhere. 

The continuum model derived from microscopic parameters has been applied to situations where also long range 
curvature effects can be significant. We have made estimates of the effects of the various spin orbit terms on the 
low energy states of fuUerenes, nanotubes, and nanotube caps. For both nanotubes and nanotube caps we find that 
All ^ IK- For nanotubes we have clarified the existent discussion and reproduced the known appearance of a gap 
for the n = states and spin-splitting for n ^ states in the energy spectrum. For nanotube caps states, we obtain 
indications that spin-orbit coupling may lead to spin-dependent emission possibilities for field-effect emission devices. 
This aspect will be investigated in the future. 

Note added. At the final stages of the writting of the present paper, two preprints [s^ l3.ll | have appeared. In 
these papers similar estimates for Ajnt ~ 10~'^meV have been obtained for the intrinsic spin-orbit coupling. Moreover 
similar discussion for the effect of an perpendicular electric field has been also discussed in|3^. Our approach is 
similar to that followed in|3^. The two studies overlap and although the model used for the a band differs somewhat, 
the results are quantitative in agreement Ag ~ lO^^meV 

The two preprints |82L Isi^ and our work agree in the estimation of the intrinsic coupling, which turns out to be 
weak at the range of temperatures of experimental interest (note, however, that we do not consider here possible 
renormalization effects of this contribution (24 156] ) . 

On the other hand, the effect of local curvature Acurv on the spin-orbit coupling has not been investigated in 
|32l l33l |. We show here that this term Acurv is as important as, or perhaps even more important than the spin-orbit 
coupling due to an electric field Ag for the typical values of £ reported. 
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APPENDIX A. TWO PARAMETER ANALYTICAL FIT TO THE SIGMA BANDS OF GRAPHENE. 



A simple approximation to the sigma bands of graphene takes only into account the positions of the 2s and 2p 
atomic levels, eg and ep, and the interaction between nearest neighbor sp^ orbitals. The three sp^ orbitals are: 

|i> ^ -L(|.) + V2b, 



(46) 




The two hopping elements considered are: 



Vi = («|Hatom|j)|i#j = 



£s - 



V2 = {i,m\Hh. 



oppingl^, f^) \m.n:nearest — neighbors — 



Vssa + 2\f2Vs'„a + 2K, 



(47) 



where i,j = 1,2,3 denote "bonding" sp^ states and n, m denote atomic sites. V\ depends on the geometry/angle 
between the bonds at each atom and V2 depends on the coordination of nearest neighbors in the lattice. Y\ and V2 
therefore determine the details of band structure for the a bands 34J. The energy associated with each "bonding" 
state (i|7iatom|*)|i=i,2,3 = (cs + 2ep) /3 is an energy constant independent of these details and not important for our 
discussion here. 

We label ai,a2,a3 the amplitudes of a Bloch state on the three orbitals at a given atom, and 
{61, 62, 63}, \h\, ^3}' {^1 1 ^2 J ^3} the amplitudes at its three nearest neighbors. These amplitudes satisfy: 



We can define two numbers, a,i = ai 



eai — Vi(a2 + as) + V26 

ea2 = ^i(ai + a3) 

ea3 = y\{a\ + 0,2) 

ebi = ^1(52 + ^3) 

eb'2 = Vi{b\+b'^) 

,bi = Yr{y[ 



b'i) 



V2b'^ 

V^b'i 
V2ai 
V2a2 
V2a3 



(48) 



a2 + 03 and bn = bi + b'2 + 63 associated to atom n. From Ea. (|48(l we obtain: 
{e-2Vi)an = V2b„ 



ie + Vi)bn = V2an + Vi 



11' \7l.—7l 



(49) 



where E„';„.-„. = (61 + 62 + ^3) 
This equation is equivalent to: 



{b'l + 62 + ^3) + (^1 + ^2 + ^3) ~ denotes nearest-neighbors. 



2Vi 



e + Vi 



n ,n.— n. 



(50) 



Hence, the amplitudes a„ satisfy an equation formally identical to the tight binding equations for a single orbital 
model with nearest neighbor hoppings in the honeycomb lattice. In momentum space, we can write: 



Vi 



(51) 



where: 



/jj = Y 3 + 2 cos(kai) + 2 cos(ka2) + 2 cos[k(ai — a2)] 



(52) 
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and ai,a2 are the unit vectors of the honeycomb lattice. 

The derivation of eauation (|51() assumes that a„ ^ 0. There are also solutions to eas. (l48|l for which a„ = at all 
sites. These solutions, and Eq.jSIJ lead to: 

ej^ = ~Vi±V2 (53) 

These equations give the six a bands used in the main text. 

In order to calculate the effects of transitions between the tt band and the a band on the spin orbit coupling, we 
also need the matrix elements of the spin orbit interaction at the points K and K' . At the K point, for instance, the 
hamiltonian for the a band is: 



( ^ 


Vi 


Vi 


V2 








\ 







Vi 














Vi 


























Vi 


Vi 












Vi 





Vi 




V 







Vi 


Vi 





/ 



(54) 



The knowledge of the eigenstates, Ea. H5U|) allows us to obtain also the eigenvalues of Ea. H54|) . The spin orbit coupling 
induces transitions from the \K,A, ]) state to the sigma bands with energies Vi ± V2 and spin down, and from the 
\K, A, [) state to the sigma bands with energies Vi /2± [%Vi) / 4, + V2 and spin up. The inverse processes are induced 
for Bloch states localized at sublattice B. 

In the hmit Vi <C V2, the a bands lie at energies ±V2, with corrections associated to Vi. The spin orbit couphng 
induces transitions to the upper and lower bands, which tend to cancel. In addition, the net effective intrinsic spin 
orbit coupling is the difference between the corrections to the up spin bands minus the those for the down spin bands. 
The final effect is that the strength of the intrinsic spin orbit coupling scales as {/^"^ /Vi){Vi/V2Y in the limit Vi ^ V2. 



APPENDIX B. MATRIX ELEMENTS OF THE SPIN ORBIT INTERACTION IN A SPHERE. 



Both the coupling induced by the curvature, eas. (|2()l21|l . and the intrinsic coupling, Ea. (|24|l . can be written, in 
a simple form, in a local basis of wavefunctions where the spin is oriented perpendicular to the graphene sheet, 
|0, 0, J-t), 1^*, J-i). Using spherical coordinates, and 0, the basis where the spins are oriented parallel to the z axis, 
I t), I I), can be written as: 



T) 
i) 



cos 



-*'^/2|0,0,±T)+cos 



-i</./2| 



(55) 



where the states |+) and |— ) are defined in terms of some fixed frame of reference. From this expression, we find in 
the basis {\A t), \A i), \B ]), \B [), basis for /C: 



s-o 



Aint cos(6') 
Aint sin(0)e*"* 



Aint sin(0)e-*^ 
-Aint cos(6') 



«AR,sin(f)cos(|) -^AR,cos2(f)e-''^ \ 



+iARsin2 (f) 



and for KJ 



i Ar sin (§) cos (f ) -i Ar sin^ (f ) e''"^ -Aint cos(6') 
V +iARC0s2(f)e*"* +zARsin(f)cos(f) -Aint sin(0)e*"^ 



iARsin(|)cos(f) 



-Aint sin(6')e 
Aint cos(6') 



(56) 



n. 



-Aintcos(0) -Aint sin(0)e-^^ -i Ar sin (f ) cos (|) Ar sin^ (f ) e"** \ 

-Aint sin(6')e*'^ Aint cos(6l) -K^Arcos^ (|) 6+*"^ +iAR sin (f ) cos (f ) 

+iARsin(|)cos(|) -zArcos^ (f)e-*^ Aint cos(0) Aint sin(0)e-''^ 

V +iARsin^(|)e*'^ -i Ar sin (|) cos (|) Aint sin(6')e*'^ -AintC0s(6l) / 



(57) 
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APPENDIX C. MATRIX ELEMENTS OF THE SPIN ORBIT INTERACTION IN A CYLINDER. 

The previous continuum analysis can be extended to nanotubes. We use cylindrical coordinates, z, cj), and, as before, 
define the spin orientations | t), | |) as parallel and antiparallel to the z axis. The matrix elements can be obtained 
in a similar way to Ea. H57|l by choosing 9 — tt/2. 



and for /C' 



n. 



n_/K, 



( Ainte-^* 

Ainte** 

-iAR,/2 -iAR/2e-"^ 

V +«AR/2e** +iAR/2 



-Ai„te** 



V +iAR/2f 



-Ainte" 




-i(p 



-iA-R/2 -iAR/2e-"^ 



,10 



+«Ar/2 
-^AR/2e+■'* 




-A 



inte 



-«Ar/2 

-^AR/2e+■'* 


Ainte^"^ 



-^AR/2e-** \ 
-«Ar/2 
-Ainte-"^ 
J 



-^AR/2e-^* 
+«Ar/2 
Ainte-*"* 




(58) 



(59) 



After integrating over the nanotube circunference J d<p the Hamiltonian above becomes: 



n_/K. 



+iAR7r \ 

-iArtt 

^^Artt 

+iAR7r 



(60) 



and for /C' 






-iARTT 

-iArtt 



-iArtt ^ 

+iAR7r 





(61) 
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